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Abstract 

The cosmic ray spectrum extends to energies above 10^*^ eV. In direct production or 
acceleration models, as well as by photo-pion interaction high energy cosmic ray flux 
must contain neutrinos and photons. The latter are absorbed by cosmic radiations 
while neutrinos are not. The need of a Neutrino Astronomy is compelling. In this 
paper a study of a detector array designed to measure horizontal r air-showers 
emerging from the ground, produced by Uj- interactions with the Earth's crust, is 
presented. Each array unit is composed of a pair of scintillator tiles mounted on a 
frame with a front field of view of about 0.1 sr, optimized to distinguish between 
up-going and down- going crossing particles by their time of flight. The detector 
array sensitivity, the size of the array and the r shower identification are discussed. 
Because of the almost complete mixing of v-t the ultrahigh energy neutrino 

tau and its minimal consequent tau-airshower rate is estimated; assuming that the 
neutrino energy spectrum follows a Fermi-like power law E^"^, the sensitivity with 
3 years of observation is estimated to be about 60 eV cm^^s~^sr~^ in the energy 
range 10^''"^'^ eV. This value would provide competitive upper limit with present 
and future experiments. We found also that, in the same time, this system can 
observe about one GZK neutrino event per km^. 
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1 Introduction 



In the rising era of UHECR astronomy, the understanding of astrophysical phe- 
nomena can shed Ught on the structure and evolution of the Universe. In par- 
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ticular, the observation of ultra-high energy cosmic rays (UHECRs) has raised 
new questions about their source and their composition. The old Ginzburg 
scheme according to which cosmic rays come from galactic sources such 
as supernova remnants or pulsars, has received new attention as they can be 
possible neutrino sources However Supernova may provide up to PeVs 

energy accelerations. Cosmic rays are mainly protons and atomic nuclei that 
are observed to strike the Earth with exceedingly large energies, but lower 
and lower fluxes. Their energy spectrum extends to energies above 10^°eV and 
is mostly well described by a power law with two kinks, a "knee" , where the 
slope steepens at energies of about 3 x lO^^eV, and an "ankle" at energies of 
3 X lO^^eV, where the spectrum flattens. It is generally accepted that most 
cosmic rays below the knee may come from inside our galaxy while highest 
energy ones, above the ankle, from extra-galactic sources. However, the origin 
of UHECR's, i.e. those above the knee, remains a puzzling mystery. Nobody 
knows with security where they are produced, and no guaranteed source has 
been detected yet. The idea that high-energy cosmic rays originating in our 
galaxy and those originating from extra-galactic sources could come from the 
same kind of source has been suggested ^ 0). The origin of UHECRs above 
the ankle could be due to particles generated in the cores of Active Galac- 
tic Nuclei (AGN), which could then produce neutrinos through interactions 
of accelerated protons and ions with an accretion disk target surrounding a 
central black hole (0). In analogy UHECR has been suspected to be produced 
by GRBs and SGRs sources. However, the UHECR flux is expected to be 
suppressed after interacting with the cosmic microwave background radiation 
when cosmic rays have energies above 4^5 10^^ eV. Such a rate drop is gener- 
ally referred to as the "GZK cutoff" (j^, following the Greisen, Zatsepin and 
Kuzmin proposal. If the GZK cutoff is not confirmed by experiments, new 
physics must be invoked to understand the origin of these extreme events, and 
the neutrinos could play a crucial role both as a source |l9[ ) as well as a probe. 



Indeed neutrinos are often associated with high energy cosmic ray fluxes and, 
since they interact only weakly, they carry directionality as a characteristic 
signature of their source. Therefore, high energy neutrino observations are 
valuable probes of new astronomy. Information is encoded in the energy spec- 
trum, arrival direction, and flavor content of the cosmic neutrinos. The mea- 
surement of the neutrino-flavor mixing rate is also important to determine the 
origin of these particles. These UHE particles may come from very rare exotic 
sources due to dark matter decay (45-) and annihilations 



Most experiments can detect neutrinos by looking for the most penetrating 
up-going lepton tracks (the muons) produced by charged current interaction 
of parental neutrinos with energies of 10^^ — 10^^ eV. Neutrino atmospheric 
secondaries from common cosmic ray flux in this energy range are investigated 
using several techniques (11; l^; 
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However a totally novel technique has been proposed based on the peculiar 
tau decay in flight: the discovery of UHE tau neutrinos at Horizons or Upgoing 
(Hortaus and Uptaus or Earth- Skimming Neutrinos) ([lil: FlQ.: .2(k .2l[ whose 
powerfuU air-shower tests the PeV-EeV tau neutrino astronomy. The up-going 
signature is noise free respect to downward ones. 

Indeed the on-going experiments, which are designed to detect particle show- 
ers with energies above 10^^ eV and are sensitive to neutrinos coming from 
an almost horizontal direction. For example, the Southern site of the Pierre 
Auger Observatory is able to detect down-going air-showers through photo- 
luminescent light at azimuthal angles of 6* < 60°, and has also a large ac- 
ceptance to the showers produced by Earth- skimming neutrinos with energies 
above lO^^eV (HIQ. 



None of the experiments devoted to the exp loration of neutrino physics or 
UHECRs, such as Amanda/Icecube Il2[ ). Antares, Magic or NuTel, are 
specifically designed to measure the flux in the energy region above lO^^eV 
;ith a la;ge duty cycle ©. The proposed detector array described here is 
designed with improved sensitivity to measure or set an upper limit on the 
flux of Earth-skimming tau neutrinos. This measurement will introduce new 
constraints on the flux of UHECRs, and distinguish between the AGN and 
the GZK production mechanisms. 

The paper is organized as following. In Section 2, the results of simulation 
studies of atmospheric hadron and inclined r showers are presented; configu- 
ration of the detector array, detector acceptance and predicted event rates are 
discussed in Sections 3-4. In Sections 5-6, the selection and the reconstruction 
of the tau showers are studied; Section 7 contains a discussion of the trigger. 
In Section 8 the results obtained with a 2-detector array prototype installed 
at the High Altitude Jungfraujoch Laboratory (3600 m a.s.L, Switzerland) are 
presented. 



2 Showers produced by UHE hadrons and z/t-'s 



The ability of detecting tau neutrinos increases significantly for neutrinos that 
enter the atmosphere at large zenith angles, at horizontal or almost below the 
horizontal direction, traversing a small section of the Earth's crust. In fact, at 
energies above 10^^ eV the entire Earth becomes opaque to neutrinos and only 
horizontal or Earth-skimming high-energy neutrinos can be detected. Neutri- 
nos coming from a vertical direction have a negligible probability to interact 
with the atmosphere (a ten meter water equivalent screen) and remain mostly 
invisible to detectors above the ground. When Earth-skimming neutrinos in- 
teract, they may generate tau whose interaction lenghts may exceed range in 
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rock of 5-10 km at energies above 10^^ eV (Fig. 1) (|2]Jl. 



The ratio of the interaction length of the tau traversing the rock divided by 
the interaction length of the neutrino, Xr-rock/^u) is about 1% in the energy- 
interval of lO^s-i^eV. 
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Fig. 1. The neutrino interaction length as a function of the energy for in water and 
rock and tau range in air, water and rock (|2lh. 



The main source of background to the neutrino flux comes from atmospheric 
protons and muons, and the most difficult experimental challenge comes from 
separating this background from inclined tau showers. The aim of the study 
presented in this section is to investigate the shower properties using a Monte 
Carlo simulation. In fact, the topology of the shower can be used to distinguish 



UHE neutrinos from the ordinary atmospheric hadron interactions ()2flt l22l:l24: 

\2 



Ordinary hadrons interact when entering the atmosphere and initiate down- 
going air-showers. At large zenith angles {9 ~ 90°) they interact at a distance 
of about 400 km from the ground, with a maximum development of the shower 
at approximately 100 km after the interaction occurs. Because of interactions 
with the atmosphere, at the ground level, most of the electromagnetic com- 
ponent of the shower disappears and only more energetic muons survive. The 
density of muons coming from protons with an initial energy of 5 x lO^^eV is 
estimated using the Aires simulation (jioh to be of the order of 10~^ muons/m^. 
On the other hand, UHE neutrinos interact weakly with protons and nuclei 
through charged or neutral currents producing leptons or neutrinos with en- 
ergies of 10^'^~^*'eV. When a Vr interacts through charged current, high energy 
r leptons are produced with long decay lengths in the range of 5-50 km which 
then originate air showers with a small transverse size (a few kms), if occurring 
near the observer. The arrival direction of the can also be inferred from the 
shape of the shower itself. 
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Simulation studies suggest that the time structure and the muon density sig- 
nificantly differs, at ground level in specific geometrical conditions, in r and 
hadron showers (see section 5). In this study only the tau hadronic decay 
modes (for instance, r — > vrvrvrz/^) have been considered, which are about 64% 
of all tau decay modes. The shower initiated by a tau lepton decaying into 
TCTnTUr has a density of about 10'^ times larger than a horizontal atmospheric 
proton shower of the same energy and zenith angle due to the long path in 
the atmosphere. 



3 Design of the detector array 



In order to detect muons produced by inclined air showers, several detectors 
pointing at the horizon and arranged in a large surface array must be as- 
sembled. The elementary modules, called towers ()27l). can be arranged in a 
grid with a spacing of distance D (Fig. 2). In a grid with detector spacing 
Z)=100 m, for example, the time correlation of the signals from the individual 
towers of about 300 ns allows the identification of the arrival direction of the 
shower. 




Fig. 2. Schematic aerial top view of a partial detector array. Four rows of twelve 
towers each are spaced by a distance D (not to scale). 

Due to the fact that the Earth-skimming neutrinos are coming from inclined 
showers, detector arrays pointing at the horizontal direction are more efficient 
than those pointing at the sky at small zenith angles. Therefore, the location of 
the detector array, pointing downward on an inclined plane or a mountain slope 
at large elevation to cover a large solid angle, can substantially increase the 
detector's acceptance and improve the rejection of background events (Fig. 3). 
Ground detector arrays have small trigger efficiency for inclined showers close 
to the 90° zenith angle. They can mostly detect horizontal r showers generated 
in the atmosphere. 

In order to measure upward moving showers coming from Uj. interactions with 
the Earth's crust, the detector must be able to distinguish between up-going 
and down-going particles, and towers can be individually used for this pur- 
pose. Each tower (Fig. 11) consists of two parallel scintillator plates placed at 
a distance of 160 cm and is instrumented with a precision timing measurement 
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Fig. 3. Schematic drawing (not to scale) of the detector array located on the 
slope of a mountain along a zenith angle of 95° (left), and pointing at the incoming 
inclined showers from Earth's skimming neutrinos (right), after traversing a distance 
L through the Earth's crust. 



device that allows to measure the time of impact of a particle passing through 
the scintillator tiles. The particle direction can be determined and up-going 
particles can be selected. The charge deposited can also be measured. Results 
of a tower prototype designed to measure the horizontal flux have been pre- 
viously reported ()27ll. Further studies to determine the horizontal flux from 
atmospheric showers have been performed at the Jungfraujoch Laboratory, 
which is located at 3600 m a.s.l., and are reported in section 8. 



The shower shape at ground level has been studied to determine the needed 
granularity of the detector array, i.e. the distance between towers. The muon 
density of a 5 x lO^^eV r shower emerging from the ground is about 10 
particles/m^, in average, and a minimum of 2 particles/m^ at the border of 
the shower; therefore, using a tile size of 20 x 20 cm^ and assuming a constant 
shower density, each tower has a probability of 40% and 8% to detect a particle 
in the core and at the border of the shower, respectively. 

Due to the relatively short longitudinal development of the shower |28|), the 
distance between the ground, from where the shower is emerging, and the de- 
tector array should be kept short and preferably about 10 km. At increasing 
ground-to-detector distances the contamination from background events com- 
ing from above the horizon also increases. No u^- regeneration from r decays 
has been taken into account (H^. The study of the shape of a r shower at 
ground level also indicates that two detector arrays located on two nearby hills 
or slopes of a mountain may improve the acceptance, i.e. using two rectangular 
arrays separated by a gap. 
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4 Detector acceptance and predicted event rate 



The flux of 10 eV is estimated to be of the order of 10 km~ yr~ sr~ in the 
GZK model and about 10^ larger in the AGN production mechanism (Fig. 4, 
left). According to these expectations, in order to detect at least 1 event per 
year, the detector acceptance, defined as detector sensitivity, must be of the 
order of 1-10 km^sr. 
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Fig. 4. Energy spectra of the predicted neutrino flux from the AGN model: MOSlB, 
P96|3l!) and 8891(0) (left); cosmological neutrino flux compared with expected 
sensitivity to tau neutrino flux for the planned projects (right). From Ref . (|33: IsT 



The probability that a lO^^eV Vr produces a 10^''eV tau emerging from the 
ground depends of the distance L traversed through the Earth's crust (Fig. 3) 
and on the tau energy loss in the rock. This probability is calculated as the 
convolution of the probability that the Vj. interacts before a distance L with 
the probability that the r reaches a distance L. Figure 5 shows, considering 
Vr flux with E^"^ power law in the energy range 10^'^ — lO^'^eV, the number 
of taus reaching the ground, which decreases as the distance L increases. For 
L=200-300 km, which corresponds to a zenith angle of 92.5°, the conversion 
efficiency varies from 1.5% to 3% as function of L along the length of the 
detector array between 10 km and 5 km, respectively (Fig. 6). The conversion 
efficiency is defined as the probability of getting a tau shower at ground level 
having a shower maximum close to the detector plane. 

The predicted event rate is calculated for a detector array located on a surface 
inclined at an angle a with respect to the ground (Fig. 9). The solid angle 
can be improved if the detector array is composed by units, each formed by 
two towers, shown in Fig. 11 and positioned nearby (i.e. 60 cm) and parallel to 
each other. In this case, if each unit has opening angles A0 = 40° and AO = 
15° and Sgff is the effective area defined like the projection of the detector 
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Fig. 5. Neutrino interactions along the Earth's crust where a tau reaches the ground 
(shaded blue). The red shaded area is the fraction of tau emerging the ground and 
decay with the maximum of its shower on the detector plane distant 10 km from 

the emerging point. For flux is considered a power law E~'^ in the energy range 
1017 _ io20ey. 
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Fig. 6. Conversion efficiency for different ground-to-detector distances (5,8,10 km) 
as function of distance L, the combined path of Vr and r through the Earth's crust. 
For v-r flux is considered a power law E~'^ in the energy range 10^'' — lO^'^e^. 
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plane to the transverse shower section, we obtain an acceptance of about 
S'e//0.15 krn^sr. Assuming that the energy spectrum of cosmological neutrinos 
follows the power law d^/dE^ = 10~^i?~^GeV~^cm~^s~^sr~^, an integrated 
flux of 3 ■ 10^km~^yr~^sr~^ is expected in the energy range 10^^ — lO^^eV, 
which corresponds to approximately 5 events per yr km^. Assuming that the 
neutrino energy spectrum follows a Fermi-like power law the sensitivity 
with 3 years of observation is estimated to be about 60 eV cm~^s~^sr~^ in the 
energy range 10^^~^° eV. This value would provide competitive upper limit 
with present and future experiments. We found also that, in the same time, 
this system can observe about one GZK neutrino event per km^. It is therefore 
possible to exclude or confirm several models ( AGN-M95 (jssl ) . MPR ljsil )) and 
determine the energy spectrum of UHE neutrinos. 



5 Identification of a tau shower 



The arrival time of the particles can provide information on the arrival di- 
rection of the shower. The ordinary horizontal hadron interaction has a plane 
wave front since it originates very far from the detector while the r shower is 
characterized by a front with a curvature radius R^IO km on detector plane 
(Fig. 7). Therefore, a time measurement of the shower development provides 
the ability to identify the direction of the reconstructed shower and the decay 
point. The muon density generated at the detector array from atmospheric 
hadron showers and from tau showers are compared using the Aires simula- 
tion (Fig. 8). In the simulation, the hadron shower arrives at a zenith angle of 
89°, while the r shower emerges from the ground at a zenith angle 6 = 95°. 
Both showers have an energy of 5 x lO^^eV and are sampled in a time window 
of 500 ns. 




Fig. 7. Schematic drawing of the time propagation through the detector array of a 
shower originating in the origin O which exhibits a curvature radius R at detector 
level; d is the distance from the axis of the shower direction. 
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Fig. 8. Muon density produced in a time window of 500 ns at the detector location 
is estimated using the Aires simulation: atmospheric (left) and r showers (right) 
have an energy of 5 x lO^^eV and are coming from the right side of the plots. 

Different strategies can be combined to distinguish a r shower from other 
particles originating from ordinary UHECRs: 

• detector array pointing at different zenith angles (i.e. 6 = 92° — 95°) and 
leaving the possibility to change the azimuth angle of each unit to enlarge 
the acceptance, (|37[ ) ; 

• measure the particle density which is different for a "new" UHECR air 
shower with a depth ^ of about 400 g/cm^ [R fa 10 km), and an "old" one 
with a depth larger than 4000 g/cm^ [R >100 km); 

• measure the direction of the shower axis and its vertex; 

• measure the ratio between the short lived electromagnetic component and 
the long-lived muons of an air shower (the possibility to distinguish between 
these two components is still under study). 



6 Tau shower reconstruction 



A study of the tau shower reconstruction has been performed using Aires (j26|) 
and Corsika (jsih Monte Carlo simulations. Tau shower Monte Carlo events 
with a zenith angle 6 = 95° have been generated and the detector array is 
assumed to lie on an inclined plane (i.e. mountain) of 45 degrees and positioned 
10 km away from the ground where the r is emerging . The conclusions are 
not strongly dependent on theta angle within 5 degrees. In the simulation we 
consider muons and electrons with energy greater than 300 MeV and 20 MeV 

^ The "shower depth" is the distance between the tau decay and the detector. 
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respectively. The identification of the shower core, if it is inside the array, is 
estimated by calculating the center of gravity of the shower with an accuracy 
of about 5 m. This is affected by an additional uncertainty of approximately 
50 m due to the direction and inclination of the incoming shower, and to the 
exponential decrease of particle density. With an angular resolution of less 
than 1°, the intersection of the shower axis with the Earth's crust can be 
measured with an uncertainty of approximately 100 m. 

The direction of the shower has been evaluated using a least squares fitting 
technique for a function describing the position {x, y, z) of the front of the 
shower and the recorded arrival time U at the detector site, with respect to 
the time of decay (to) of the tau, by minimizing 

= ^^=i{xil + Uim + Zik - c(ti - to))^ 

where the sum includes all the towers and I, m, k are the direction cosines. 
Using the muons in the shower, the resulting zenith resolution is estimated to 
be of the order of 0.5°. This result is independent on the size of the subsection 
of the detector array, when an array grid of at least 100 x 100 m^ with spacing 
D =25 m is set. It should also be remarked that each tower in the detector 
array is "blind" to particles coming from a direction that deviates of more 
than 7.5° from the pointing axis. Furthermore, if the shower is contained in 
the detector array, the shower front can be used to estimate its radius, i.e. 
distance of the decay point of the tau (or shower depth), and thus gain a 
further rejection of background events. 

The possibility of detecting showers with the core outside the array will be 
considered due to the large acceptance of units, although the contamination 
due to background events may be larger. If the core is outside of the array 
the reconstruction efficiency is estimated to be of the same order even if not 
accurately evaluated. 




Fig. 9. Schematic drawing of the particle shower front impinging on the detector 
array. 
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7 Trigger 



The detector array is composed by several units spread over an area of 1km x 1km. 
The remote location and the lack of easy communication links requires reliable 
low power electronics powered by solar panels, and trigger logic at each tower. 
Trigger algorithms, operating within each tower and between different towers, 
must be developed to suppress lower energy cosmic ray showers and retain the 
events of interest. A hierarchical event trigger can be used in order to keep 
the event rate within the limits of the data acquisition system. Correlation 
of more towers can provide a further rejection of air showers. Furthermore, 
the particle multiplicity from UHECRs is larger than in atmospheric showers. 
The timing information of the front of the shower can also be used to discrim- 
inate between the events of interest with short radius {R ~ 10km), and the 
atmospheric air showers with much larger radius {R > 100km). 

Following these criteria it is possible to address different physics issues and 
select r showers, muon bundles or large zenith angle atmospheric interactions. 



8 Test results with a detector prototype at Jungfraujoch Station, 
Switzerland 



Measurements were performed at the High Altitude Jungfraujoch Station 
(3600 m a.s.l.) to understand detector characteristics and performance. These 
studies are also aimed at the understanding of the possible sources of back- 
ground from inclined atmospheric showers at large zenith angles. The envi- 
ronment of this test was intentionally chosen to be worse than at the final 
experimental proposed site. At an altitude of 3600 m a.s.l. the electromag- 
netic component of vertical air showers is larger than at sea level, increasing 
the probability of contamination from background events. Studies of the de- 
tector performance with one module called tower at sea level were reported 



earlier (1271). 



Towers with different tile sizes have been installed: two are instrumented with 
tiles of dimensions 12.5 x 12.5 x 2 cm^ and are placed parallel to each other 
about 50 cm apart, while another tower has tiles of 20 x 20 x 1.4 cm^ and 
was installed at a distance of 20 m from the other two (Fig. 10). The use of 
tiles with different sizes is aimed at the optimization of detector acceptance 
and time resolution. NIM and Camac modules with low threshold discrimina- 
tors, high resolution Time-to-Digital converter (TDC), and Analog-to-Digital 
converter (ADC) were used to trigger the tracks and measure their time of 
flight. The readout electronics used in these measurements will be replaced by 
a customized board installed on each tower for the final design of the detector 
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array; using TDC-GPX by Acam, with its 80ps resolution, and MATACQ by 
DAPNIA/CEA, capable of a continuous sampling at IGHz for 4/8 channels, 
for a memory depth of about 2500 points, it will be possible to implement 
a fully multi-hit solution for time, charge and pulse shape analysis within a 
2.5fis window. 

Each tower is composed of two parallel tiles (CI and C2) of Kuraray scintilla- 
tor mounted on an aluminum frame (Fig. 11), and are placed at a distance L 
which can be varied from 1.4 cm to 200 cm. The distance is set at L=160 cm 
to distinguish between up-going and down-going tracks and to reduce the 
background generated by "mini-showers" , due to particles traversing the scin- 
tillator near the edges (j27|). The direction of the particles traveling through the 
tower can be determined by the measurement of the time of flight between the 
two tiles. The scintillators are read by a Hamamatsu H5873 photomultiplier 
tube (PMT) directly coupled to the scintillator through an optical connector, 
and without a light guide in order to minimize the time dispersion of the sig- 
nal. The PMT has a low-voltage power supply and a transit time spread of 
less than 0.8 ns. 




Fig. 10. Two towers located at the Jungfrajoch High Altitude Station. 



8. 1 Detector performance 



A measurement aimed at determining the time resolution was performed using 
a single tower pointing at a zenith angle 6* = 0° and with a distance between 
tiles L =160 cm; the difference of the signal arrival time at each of the two 
tiles, Ati2, was recorded (Fig. 12). The intrinsic time resolution of the tower is 
measured by the spread of the distribution a ~1.2 ns. A similar time resolution 
was obtained when the two tiles were placed at a distance L =1.4 cm. 

In order to isolate the particles arriving longitudinally along the axis of the 
tower, selection cuts are applied to the charge collected by the two photomul- 
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Fig. 11. Schematic view of a tower (not to scale). The up-going, down-going, and 
vertical tracks shown can be distinguished using the time-of-flight information. 
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Fig. 12. Difference of time-of-flight (Ati2) between CI and C2, when the 20 x 20 cm^ 
tiles are 160 cm apart. 

tipliers. Orthogonally-traversing tracks deposit in the scintillator tiles approx- 
imately the same amount of energy, within Landau fluctuations. Asymmetri- 
cally deposited charges characterize background tracks: very low or very large 
amounts of deposited charge are associated to tracks traversing near the edges 
of the scintillator or to multiple tracks, respectively. These cuts were deflned 
using the charge of vertical equivalent muons (CVEM). Data were collected 
using a coincidence of the signals in CI and C2 with a gate of 70 ns. The light 
emitted by the scintillator has a cylindrical symmetry and the time resolution 
is determined by the flrst photons that reach the PMT, without any reflection 
on the edges. Photons hitting the corners of the tile will arrive at the PMT 
at a later time and, due to the reflections, will have a smaller charge at the 
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PMT's output. 



8.2 Measurements 

A measurement of the horizontal cosmic ray flux was performed using one 
tower with the tiles placed at a distance L=160 cm, and pointing at a zenith 
angle 95°. Thanks to its good time resolution, each tower is capable of dis- 
tinguishing the particle direction between up-going and down-going tracks by 
measuring the time- of- flight between the two tiles. The two peaks at -5 ns and 
-|-5 ns correspond to up-going and down-going tracks, respectively (Fig. 13). 
The contamination of the events outside the two peaks is due to parallel tracks, 
most hkely vertical, or to tracks where the particle hits an "inefficient region" , 
i.e. a corner, of one of the tiles. 

The contamination in the flux above 90 degrees is due to muon decays and 
can be eliminated by absorber material (lead) or, alternatively, by correlating 
information from distant towers. The suppression provided by the lead block 
is different for up-going and down-going tracks, and it is due to the shadow of 
the mountain which remove part of the down-going track flux. 

A measurement of the cosmic ray flux was performed using two towers pointing 
at different zenith angles and compared to results from other experiments 
(Fig. 14). The detector sensitivity is 10~^cm~^sec~^sr"^ at ^ = 98°. 



9 Conclusions 

The study and the implementation of a detector array designed to measure 
the flux from UHE neutrinos has been discussed. Results using a detector pro- 
totype have also been presented. In the hypothesis of full mixing of z/^ — > z/^, 
after two years of data-taking with the full detector array, about 10 events 
could be collected in the energy range 10^^ — 10^°eV from AGN tau neutrinos, 
or an upper limit on the cosmic ray flux, as well as measurements of muon 
bundles at large zenith angles, could be achieved for other scenarios. Assum- 
ing that the neutrino energy spectrum follows a Fermi-like power law E^"^, 
the sensitivity with 3 years of observation is estimated to be about 60 eV 
cm^^s^^sr"^ in the energy range 10^^~^° eV. This value would provide com- 
petitive upper limit with present and future experiments. We found also that, 
in the same time, this system can observe about one GZK neutrino event per 
km^. Even if we have not investigated accurately cases in which the core of the 
shower is outside of the array, we expect the reconstruction efficiency not to 
decrease signiflcantly. The proposed detector array is complementary to other 
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Fig. 13. Time-of-flight difference between two tiles Ati2 in a tower pointing at 
a zenith angle 95° when a lead block of 3-cm thickness is placed in front of one 
tile (yellow shade) and no lead is present (red shade). The two peaks at -5 ns and 
+5 ns correspond to up-going and down-going tracks, respectively. The peaks at -3, 
and 3 ns are due to parallel tracks, most likely vertical, where one of them hits 
an "inefficient region", i.e. a corner, of one of the tiles. 

experiments in this energy range, like Amanda or Antares. It is also more 
efficient than those experiments using "double-bang" or "lollipop" strategies 
which suffer from the inefficiency due to the need of detecting the hadronic 
shower from the initial z/,-, where the r track and a second cascade must be 
close to the detector. In terms of energy range, it is comparable to the Auger 
and NuTel experiments but the former has no large detection efficiency for 
taus emerging from the ground and the latter has only 10% of the duty cycle. 
A recent proposal that uses a fluorescence technique (0) will be competitive 
with this project in terms of detected events if two new telescopes pointing at 
a long chain of mountains are built. 



Acknowledgments 



During the writing of this paper we have beneffied from discussions with 
several colleagues; among them J. Russ, J.F. Beacom, M.A. Huang . 



16 



Cosmic ray flux 



JO 

L 

'e 

> 

"5 -4 
=10 

0) 



£10 

3 

lio' 
< 



10 



10 



10" 



















■ 

1 


1 














• 






• This Test at 3600 m 
A Nevod, p^> 7 GeV/c 

A Karmakar et al. 
■ Wilson et al. 

Jakeman et al. 
□ Flint Nath 
A Crookes Rastin 




c 

A 
























A 

A 

i 






, , , 


, , , 


, , , 


A 

1 1 il 





20 



40 



60 



80 



100 
(degrees) 



Fig. 14. Measurements of fluxes at several zenith angles in the High Altitude 
Jun gfra ujqch Station compared with results from other experiments at sea level, 



References 



[I] V.L. Ginzburg and S.I. Syrovatskii, "Origin of cosmic rays", The Macmil- 
lan Co., New York (1964). 

[2] S. Nagataki, Astrophys. J. 600 (2004) 883. 

[3] W. Bednarek, Astron. Astrophys. 407 (2003) 1. 

[4] A.D. Erlykin, T. Wibig and A.W. Wolfendale, New Journal of Physics 3 
(2001) 18. 

[5] J. Arons, Astrophys. J. 589 (2003) 871. 

[6] W. Bednarek and R. J. Protheroe, Mon. Not. Roy. Astron. Soc. 302 (1999) 

373, astro-ph/9802288, 
[7] V.K. Dubrovich, D. Fargion, M.Y. Khlopov. (Astropart. Phys. 22 (2004) 

183. 

[8] K. Greisen, Phys. Rev. Lett. 16 (1966) 748; G.T. Zatsepin and V.A. 

Kuzmin, JETP Lett. 4 (1966) 78. 
[9] D. Elsaesser and K. Mannheim, Astropart. Phys. 22 (2004) 65, 

|a'stro- ph/0405347, 

[10] M. Tekedat et al. (Agasa Collaboration), Astropart. Phys. 19 (2003) 447. 

[II] P. Niessen (for the Amanda Collaboration), |astro-ph/0306209j 
[12] For example, http://icecube.wisc.edu 

[13] T. Abu-Zayyad et al. (High Resolution Fly's Eye Collaboration), 



17 



|astro- ph/0208243, 

[14] The Pierre Auger Project Design Report by Auger Collaboration, 

FERMILAB-PUB-96-024. 
[15] D. DeMuth et al., Astropart. Phys. 20 (2004) 533. 
[16] X. Bertou et al., Astropart. Phys. 17 (2002) 183. 
[17] M. M. Guzzo and C. A. Moura, hep-ph/0504270, 

[18] Fargion, A. Aiello, R. Conversano, 26th ICRC, Salt Lake City,USA, 

HE.6.1.10, pp. 396398(1999) 
[1 9] D. Fa rgion, B. Mele, A. Salis. Astrophys. J. 517 (1999) 725, 

|astro-ph/ 9710029t 

[20] D. Fargion, Astrophys. J. 570 (2002) 909. ^ 

[21] D. Fargion, et. al; Astrophys. J. 613 (2004) 1285, |hep-ph/0305128 



[22] J.L. Feng et al., Phys. Rev. Lett. 88 (2002) 161102. 
[23] E. Zas, astro-ph/0504610, 

[24] J.F. Beacom et al., Phys. Rev. D 68 (2003) 093005.^ 

[25] Z. Cao, M.A. Huang, P. Sokolsky and Y. Hu, |astro-ph/0411677| 
[26] S. J. Sciutto, Aires, "A system for air shower simulations" (v 2.8.0) March 
2005. 

[27] M. lori, A. Sergi and D. Fargion, |as'tro-ph/04 09159. 
[28] T. Gaisser and AM Hillas, Proceedings of the 15th ICRC, Plovdiv, Bul- 
garia, 8 (1977) 353. 
[29] J.F. Beacom, P. Grotty and E. W. Kolb, Phys. Rev. D 66 (2002) 021302. 
[30] K. Mannheim, Astropart. Phys. 3 (1995) 295. 
[31] R. J. Protheroe astro-ph/9607165 

[32] F. W.Stecker and H. Salamon, Space Sci. Rev. 75 (1996) 341. 

[33] Dmitry, V. Semikoz and Gunter Sigl, JCAP 0404 (2004). 

[34] F. W. Stecker, C. Done, M. H. Salamon, and P. Sommers, Phys. Rev. 

Lett. 66 (1991) 2697. 
[35] K. Mannheim, Astropart. Phys. 3 (1995) 295. 

[36] J. P. Rachen, R. J. Protheroe, and K. Mannheim, Phys. Rev. D 63 (2001) 
023003. 

[37] D. Fargion et. al., "astro-ph/0501079" 

[38] D. Heck and J. Knapp, "Extensive Air Shower Simulation with Corsika", 

V 6.204, July 2005. 
[39] Nevod Collaboration 29th Int Conference Pune 2005 101-104. 
[40] N.L. Karmakar A. Paul and N. Chaudhuri, Nuovo Cimento B17 (1973) 

173. 

[41] B.C. Wilson, Can. J. Phys. 37 (1959) 19. 

[42] D. Jakeman et al. Can J. Phys. 34 (1956) 432. 

[43] R.W. Flint, R.B. Hicks and S. Standil, Can J. Phys. 50 (1972) 843. 

[44] J.N. Crookes and B.C. Rastin, Nucl. Phys. B39 (1972) 493. 

[45] P. Bhattacharjee, G. Sigl, Physics Report 327(2000) 109-247 



18 



